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The Dirac spectrum of QCD with dynamical fermions at nonzero chemical potential is 
characterized by three regions, a region with a constant eigenvalue density, a region where 
the eigenvalue density shows oscillations that grow exponentially with the volume and the 
remainder of the complex plane where the eigenvalue density is zero. In this paper we 
derive the phase diagram of the Dirac spectrum from a chiral Lagrangian. We show that 
the constant eigenvalue density corresponds to a pion condensed phase while the strongly 
oscillating region is given by a kaon condensed phase. The normal phase with nonzero chiral 
condensate but vanishing Bose condensates coincides with the region of the complex plane 
where there are no eigenvalues. 
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I. INTRODUCTION 

Because of the Banks-Casher relation jl| between the chiral condensate and the spectral density 
of the QCD Dirac operator, lattice QCD studies have had a long lasting interest in the Dirac 
spectrum. The eigenvalue correlations have been well understood for an (anti-)Hermitian Dirac 
operator. In the chirally broken phase, the eigenvalue density around zero is nonzero in the 
thermodynamic limit resulting in a discontinuity of the chiral condensate. In the chirally restored 
phase, a gap develops around zero so that the condensate remains a smooth function of the mass 
in the thermodynamic limit and therefore vanishes in the massless limit. 

At nonzero quark chemical potentialjU^the Dirac spectrum also received considerable attention 
in the lattice QCD literature , [lOj. In quenched lattice simulations it was 



found that the Dirac eigenvalues are distributed more or less homogeneously inside a strip that is 
symmetric about the imaginary axis. In particular, this implies that the quenched chiral condensate 
vanishes in the massless limit for arbitrarily small chemical potential Q|. On the other hand, 
because at very low temperatures the free energy of QCD with dynamical quarks does not depend 
on the quark chemical potential below a third of the nucleon mass, chiral symmetry must remain 
broken at small nonzero chemical potential if it is broken at zero chemical potential. Naively, 
the fermion determinant will make matters worse (an even smoother behavior of the condensate 
at m = 0) because eigenvalues will be repelled from the position of the mass js], but as will be 
discussed below, it is precisely the phase of the fermion determinant that will lead to a nonzero 
chiral condensate. 

A partial resolution of this paradox appeared soon after the publication of 2]. Based on a one- 
dimensional U{1) lattice model it was concluded that phase fluctuations of the quark determinant 
are essential to obtain a nonzero chiral condensate. In the same paper, it was found that the width 
of the strip of eigenvalues is equal to the quark mass when fi = This has been verified by 

lattice simulations BQ. Recently, this „snlt was derived fen a cM.al La,..a.„ian fa the Di.ac 
spectrum of quenched QCD However, the paradox of how a /u-dependent spectral density 

can lead to both a /^-independent partition function and a discontinuity of the chiral condensate 
remained unsolved and was later coined [3] as the "Silver Blaze Problem." 

The resolution of the "Silver Blaze Problem" became possible because of two further develop- 
ments. First, the idea to explain the failure of the quenched approximation at /i 7^ in a random 
matrix model with the symmetries of QCD at /i = [l^. Second, the introduction of a mathe- 
matically simpler version of this model that made it possible to derive exact analytical expressions 



for the spectral density of the 



results it was shown in 



16|, 



QCD Dirac operator with dynamical quarks Il5l ]. Using these 



Q,Q]. 



13] that a strongly oscillating contribution to the spectral density is 
responsible for the discontinuity in the chiral condensate. 

The random matrix results are valid for QCD in the microscopic limit [3] where the 
combinations 

m = mSy and fi"^ = fi'^F^V (1) 

are kept fixed in the thermodynamic limit (Here, m denotes any of the quark masses or the mass 
scale for which the Dirac spectrum is calculated, fi is the chemical potential, T, is the magnitude 
of the chiral condensate and F is the pion decay constant). Since in this limit, the contribution of 
the zero momentum modes factorizes from the partition function [2^, it is also equivalent to the 



e regime of chiral perturbation theory 
density 



domain 
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21I ]. Random matrix predictions for the quenched spectral 



have been compared successfully to lattice gauge simulations for the microscopic 



26|. 



In this paper we show that the microscopic results for the unquenched Dirac spectrum are 
characterized by three different regions or "phases": a phase where there are no eigenvalues, a 
phase where the eigenvalue density is equal to the quenched result, and a phase where the Dirac 
spectrum is strongly oscillating. The aim of the present paper is to derive the phases of the Dirac 
spectrum from a mean field analysis of the chiral Lagrangian for the Dirac spectrum of QCD at 

In the mean field limit one evaluates the extrema of the static part of the chiral Lagrangian so 
that the leading order contribution in the p-expansion of chiral perturbation theory is exactly given 
by the large m and large ft limit of the microscopic limit of QCD. The classification of the spectral 
density of the Dirac operator is therefore valid beyond the microscopic regime up to energies that 
are well below the QCD scale. Alternatively, the locus of these phases in the complex eigenvalue 
plane can be obtained [l^ from the asymptotic properties of the complex Laguerre polynomials 
used to solve the corresponding random matrix model. 

In section |TI] we review the calculation of the phase diagram of the quenched Dirac spectrum. 
The phase diagram of the Dirac spectrum for QCD with one dynamical flavor is discussed in section 
IIIH and the dependence of the free energy on the number of replicas is discussed in section llVi 
Concluding remarks are made in section IIV[ A brief discussion of some of the results of this paper 
appeared in the proceedings [2^. 
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II. PHASE DIAGRAM FOR THE QUENCHED DIRAC SPECTRUM 



In this section we review the mean field analysis of the phases of the quenched Dirac spectrum 
originally given in 



111 ] using the supersymmetric method. 



spectral density by means of the replica trick given by 



1 1 
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n this paper we evaluate the quenched 



28| 



/,^/=°(z,z*;^) = lim--9,9,.logZ"(z,z*;M), 
n^o V n 

Z''{z,z*;n) = (det"(L' + 7o/i + z)det"(Z)-7o// + z*)), 



(2) 



where (. . .) denotes the quenched ensemble average. In the mean field limit the result for 
p^f^^(^z, z*; i^i) before taking the replica limit, n —)■ 0, does not depend on n and the limit n ^ 
can be taken trivially. Therefore it is sufficient to analyze the partition function for n = I which, 
as is illustrated by the above rewriting of the partition function, is equivalent to QCD at nonzero 
isospin chemical potential with complex quark masses. The replica argument can be made rigor- 
ous in the microscopic limit through the Toda lattice equation which links the generating 
functionals, with different replica index n. Using the Toda lattice equation one can derive an 
exact identity 



15 



22] for the quenched microscopic spectral density in the sector with topological 



charge ly 



P 



''f=^{z,z*;fi) = hz*Z^=\z,z*;fi)Z^=-\z,z*;fi). 



(3) 



Since the microscopic limit of the partition function with a pair of conjugate bosonic quarks, 
Z"=-^(z,2;*;/i), given by 



15 



21 



7n=—l 



iz,z*;(l) = \z 



2u 



1 



■ exp 



(z'^ + z*^) 

1^2 



(4) 



is a smooth function of z, the phase structure of the spectral density is given by Z'^^^{z, z*; fj,) (in 
agreement with the conclusion from the mean field replica trick). 



We will examine the partition function {z,z*;fx) to leading order in chiral perturbation 



theory. The corresponding chiral Lagrangian is given by 



11 



3ii 



(5) 



with U £ U{2) and covariant derivative defined by 



V^U = d^U -i[B^,U] 



(6) 
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FIG. 1: Phase diagram of the quenched Dkac spectrum in the complex z plane. The support of the 
quenched spectrum is between the two solid vertical black lines. The labels N and tt refer to the normal and 
pion condensed phase, respectively. 



with 




By = (5^,0 and M = . (7) 




We can neglect higher order terms in /i as long as <^ F. In the mean field limit, the phase 
structure of the partition function is determined by the static part of the Lagrangian given by 

C,tat = ^F^Tt[U,B][U\B] - ^STVM(C/ + C/t). (8) 

Notice that l^/Zstat only depends on the microscopic variables In order for the mean field limit 
to be valid for the calculation of bulk quantities, we also require that m, \z\ <C F'^/T,. 
The ground state of this Lagrangian is easily found by making the ansatz 

(cos a sin a \ 
(9) 
— sin a cos a I 

This results in the effective Lagrangian 

£3^^^ = -2n^F^ sin^ Q - S(z + z*) cos a. (10) 
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The minima of this Lagrangian are given by 



sina = and cosa = — . „ ^„ — (11) 

resulting in the free energies 

£1, = + z*) and LI,,, = -2^?F^ - ^^^^^tpP\ (12) 

respectively. The boundary between the two phases is given the condition 

S(z + z*)=Vf2. (13) 

This condition can easily be understood physically. Since the mass of the charged Goldstone 
bosons is given by 

rnl ^ ii±|^. (14) 

pion condensation occurs when ^2 > (m^/2)2 = (z + z*)S/(4F2). 

A schematic phase diagram of the quenched Dirac spectrum is shown in Fig. [H The normal 
phase (N), where sina = 0, and the pion condensed phase, where cosq = S(z + z*)/(4/i^F^), occur 
for Kz + z*)S/(4F^)| > /i^ and \{z + z*)Ti/ {AF'^)\ < /i^, respectively. Inside the pion condensation 
region the eigenvalue density (per unit area) is constant 

1 

Pn^~^{z,z*;fi) = --d,d^*Cl,,, = j^j^, (15) 
while outside this region the eigenvalue density is zero p{\x\ > 2yU^F^/Il) = —dzdz*jC,f,,, = 0. 



III. PHASE DIAGRAM OF THE DIRAC SPECTRUM FOR QCD WITH DYNAMICAL 

QUARKS 

In this section we discuss the phase diagram of the Dirac spectrum of QCD with dynami- 
cal quarks at nonzero chemical potential. We will do this by relating our results to strangeness 
condensation and pion condensation at nonzero quark chemical potentials, a connection that was 
mentioned in [l^ and discussed briefly in the proceedings 



27| 



The discontinuity of the chiral condensate arises due to the contribution from the strongly 
oscillating terms in the spectral density Q, [l^]. These terms become exponentially large with 
increasing volume, but are exponentially suppressed elsewhere {l4,[l^. Here we will show that this 
region can be identified from the mean field limit of the chiral Lagrangian and it therefore shows 
that our picture of the Dirac spectrum is valid beyond the microscopic domain. 
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As in the quenched case we use the rephca trick to express the unquenched eigenvalue density 

p^/(z,z*,m;^) = hmi-5,a,.logZ"'^/(z,z*,m;/i), (16) 
n^O V n 

in terms of generating functionals for the eigenvalue density 

Z'''^f{z,z*,m;n) = {dee{D + jofi + z)dee{D--fofi + z*)det^f{D + jofx + m)). (17) 

Below we will see that in the mean field limit of the chiral Lagrangian, the free energy does not 
always scale with n as in the case of the phase quenched partition function. Therefore, the naive 
replica trick does not necessarily work in this case. However, it turns out (see section lIVp that 
the boundary of the the oscillating region as well as the boundary of the region with a constant 
eigenvalues density do not depend on n. Therefore, these boundaries can be obtained from a 
calculation for n = 1 which will be discussed in the remainder of this section. In the microscopic 
limit this also follows from an exact identity (that can be derived from the Toda lattice equation) 
between the spectral density and the ratio of the the following partition functions (recall that the 



microscopic variables defined in ([T|) are denoted by a hat) [15 1 



p^f{z,z\rn;fs) = '^{m' - zY^ Z^^Hz, z*; A) ^""""^N^j^''^ ■ ^ 

Here, the label Nf denotes the number of fermionic quark flavors with mass m and the value of 
the label n stands for the number of quark pairs with mass z and z* (a negative value of n should 
be interpreted as n pairs of bosonic quarks). Since both the microscopic limit of Z'^^~^{z, z*; fi) 
given in and the Nj flavor partition Z^f {m; p,) = Z^i{rn) do not show a phase transition in 
the chemical potential-mass plane, the phase structure of the spectral density is determined by 
Z^f'"'^^{rh,z,z*;fi). Note that the oscillating nature must come from Z^f''^^^{m,z,z*;fl) since 
both Z"'^~^ {z , z* ; ft) and Z^f{m) are real and positive. The phase structure can be analyzed in 
the thermodynamic limit which is smoothly connected to the strong non-Hermiticity limit (/i S> 1) 
of the microscopic domain. In the large /t limit with rn ~ /l^ and |z| ~ /i^ we have that 

p^/(z,r,m;/i)e2^'+^''/(2A^)+^/'^ ~ Z^/'"=i(m, z, r ; /x). (19) 

Below we will evaluate the partition function Z^/'"^^(m, z, z*; //) in the mean field limit. The 
phase structure does not depend on the number of flavors, and we will only study the case Nf = 1. 
Alternatively, one could analyze the strong non-Hermiticity limit of the spectral density [l^ where 
the phases are determined by a competition between the quenched and unquenched contribution, 
and whether a saddle point of the expression hits the boundary on the integration interval. 
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A. Low Energy Limit of QCD Partition Function 

The QCD generating functional with n = 1 for the Nf = 1 spectral density is given by the 
partition function 

Z^-f=^'"=^(m,z,z*;/i) = {det{D + fijo + m) det{D + fi-fo + z) det{D^ + ^xJo + z*)) 

= (det(L' + /x7o + m)det(L> + ^7o + z)det(D-/i7o + 2;*)). (20) 

We can interpret the flavor with mass m as the strange flavor and z and z* as the up and down 
flavors. The isospin and strangeness chemical potential are then given by 

HI = 2n, fis = f^- (21) 

The quark chemical potential /x of the generating functional for the unquenched eigenvalue density 
thus couples to the charges carried by the pions and the kaons. That is why the chiral Lagrangian 
can give us nontrivial information about the Dirac spectrum even though the Goldstone bosons 
have zero baryon charge. The masses of the the Goldstone bosons are equal to 

m^ = ReS(m + z*)/F2, ml = ^{z + z*)/F^ . (22) 

For m — i- CX3 the partition function (1201) reduces to the quenched case discussed in the previous 



section. In this case it was found 
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13( 1 that the phase with the pion condensate corresponds to the 
region where the eigenvalue density is nonzero whereas the region with vanishing eigenvalue density 
is in the normal phase (see eq. (llSp ). As we lower the strange quark mass m, at some point the 



mass of the kaon is less than that of the pion, which leads to the formation of a kaon condensate 
Given m, we will see below that the values of z and z* for which the kaon condensate is nonzero will 
correspond to a region in the complex eigenvalue plane where the unquenched eigenvalue density 
is oscillating with and amplitude that grows exponentially with the volume. 

The low-energy limit of the partition function (j20p for real z was analyzed in [34] . Here we 



extend the arguments to the complex plane. The chiral Lagrangian for the Goldstone fields is 
given by the same expression as in the phase quenched case, 

Ces = ^F^TrV^t/V^C/t - ^STrM(C/ + ^7^), (23) 

with covariant derivative defined by 

V^U = d^U -i[B^,U] (24) 
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but in this case 

Bu = di,^odmg{^fiB + ^fJ'i,^fJ'B - ^fJ'i,^PB + f^s), (25) 
M = diag{z, z*,m), (26) 

and U € SU{3). The phases are determined by the vacuum states that minimize the free energy. 
In the mean field limit it is sufficient to consider the zero momentum modes or constant fields, and 
we only have to minimize the static part of the effective Lagrangian given by 

>Cstat = ^F^Tr[U,B][U\B] - isTVM(C/ + C/t). (27) 

B. Vacuum Configurations 

The ground states of (|27|) can be parameterized by the ansatz 

U = Ry{-/3)R,ia)RyiP), (28) 

where Rx{P) is a rotation by /? about the x-axis, etc.. Because of (12 1|) we have Bq = fi diag(l, —1, 1) 
and the static part of the Lagrangian becomes 

^stat = — 2;U^F^(1 — cos^ a) — T, [{z — m) cos^ /3(cos a — 1) + z + z* cos a + mcos a] . (29) 

The ground state is given by the minimum of the real and imaginary parts of the Lagrangian. 
For definiteness we take m > and x = Re(z) > 0. For a = the value of (i is not fixed. This 
phase, where U = 1, is the normal phase with free energy given by 

/:lt = -S(m + z + z*). (30) 

For a 7^ 0, the minimum of >Cstat is either at /? = or at /3 = 7r/2. For /? = we find 

-^« = 4;^- 
This is the pion condensed phase with free energy given by 

£^ = -2n^F'^{l - cos^ a) - T.{z + z*) cos a - Em = -Sm - 2^1'^ ^ . (32) 

This expression also enters in the relation between the spectral density and the partition function 
in the strong non-Hermiticity limit (see Eq. (jl9p ). and can be concisely written as 

p^=\i,r,m;A)e-^^^'-* ~ Z^/=^'"=i(m, i, r ; A). (33) 
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FIG. 2: Phase diagram of the Dirac spectrum in the complex z plane for one dynamical flavor. The 
support of the quenched spectrum is between the two vertical solid black lines. Unquenching introduces the 
oscillations within the ellipses which intersect the z-axis at x = ±m and x — ±(|F^^^/S — m) (dashed 
curve). The labels N. tt, K refer to the normal, pion and kaon condensed phases of the generating functional 
for the eigenvalue density. In the region between the triangle and the ellipse the spectral density oscillates 
with an amplitude that is exponentially suppressed with the volume. 

We introduced >C^tat ^ formal expression for £stat without taking into account the constraint 
I cosaj < 1. Since Z^/^^'"^^(m, £, z*\ fi) ~ exp[— in the pion condensed phase, we observe 
that the eigenvalue density in this region is of order 1. Alternatively, this relation is necessary 
to obtain an eigenvalue density that does not depend exponentially on the volume the in pion 
condensed phase. 

The mean field result for the eigenvalue density (per unit volume) is thus given by 

For P = Tr/2 the static free energy is minimized by 

cos a = with Re S(m + z*)/{4fi^F^) < 1. (35) 
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This is the kaon condensed phase with free energy given by 

= -2^2f2(1 - cos^ a) - T.{m + z*) cos a - Sz = -Sz - 2fi^F^ + ^ g^2^2 ^ • 

(36) 

The condition in (j35|) arises from the requirement that the real part of the saddle point is inside 
the integration domain. 

The phase boundaries occur where the real part of the Lagrangian for two different phases 
coincides. Notice that x > 0. The boundary between the normal phase and the pion condensed 
phase is at 

(z + z*)S = Vf^ (37) 
The boundary between the pion condensed phase and the kaon condensed phase is given by 



"^^ + g^2p2 = + Re 



(z* +m)2s2 
8/i2F2 



(38) 



With z = X + iy this can be written as 



23;S + = -(x + m)S H „ „ — , (39) 



in agreement with [3]. The boundary between the normal phase and the kaon condensed phase is 
given by 

- 2,?F^ - (x + m)E + ^''t'^tf = + - Vi^')'. (40) 



8/i2F2 ^ ' J ' g^2^2 8^2^2 

This implies that kaon condensed phase dominates in the region 

2/i^F2 < < - mS, < 4/^^^^ - (x + m)S 

2 2 y^S2 , ^ (a; + m)2s2 x'^T? , , 

mS < xS < 2^2^2^ < -(x + m)S + ^ g^,^, + 2xS - (41) 

and in a similar region that is given by the mirror image of this region with respect to the y-axis. 
Because of the relation (1331) the spectral density in the kaon condensed phase is given by 

p^/=^(z, z*, m; A) ~ 6-^^^"?-*+^^='-'. (42) 

This takes into account the requirement that | cosa| < 1 in the condensed phase. The expression 
Re[£^at] < R'6[£^tat] is a curve that splits the kaon condensed region into two parts. Inside 
the ellipse, the eigenvalue density grows exponentially with the volume. Since the free energy 



12 



in the kaon condensed phase is complex, cf. (j36p . the eigenvalue density oscillates with a period 
proportional to the inverse volume. Outside of the ellipse defined by the second equation in ()4ip 
the oscillations have an exponentially small amplitude with increasing volume. 

Below we will argue that the boundaries between the different phases of the Dirac spectrum do 
not depend on the number of replicas. However, for n = 1 the mean field result for the free energy 
given by (j36p results in a vanishing eigenvalue density in the oscillatory region. This implies that 
the resolvent for n = 1 is discontinuous across the phase boundary between the pion condensed 
phase and the kaon condensed phase. We have shown numerically that the zeros of the n = 1 
partition function are located on the same ellipse. 

The phase diagram in the xy-plane is shown in Fig. [2j The kaon condensed region is split into 
two parts by the dotted curve which intersects the x-axis at xT, = - imS. Only the region 

that is inside the ellipse shows exponentially large oscillations. The three phases intersect each 
other at 

xH = 2/^2^2, = 2fi^F^ - mS. (43) 



For y = we can distinguish the following phases 321 ]: 

• a normal phase for 2/i^F^ < xT, < mS and 2mT, < Afi'^F'^ < (x + m)S. 

• a pion condensed phase for xT, < 2//^F^ < mS and xH < mS with xS > 2//^F^. 

• a kaon condensed phase for Afj^F'^ > (x + m)S > 2mS. 

As we have seen in Fig. [21 the kaon condensed phase is split into two parts by the the curve 

^ 8 9 „9 mT, , . . 

xE = -^2^2 _ _^ (44) 

with exponentially large (small) oscillations in the region to the right (left) of this curve. Figure 
[3] shows a phase diagram of the Dirac spectrum in the fi^ — x plane for y = and m = 5. For 
fjpF"^ < mS/2 kaon condensed phase is absent. 



The border lines between the different phases can be easily understood physically. Pion conden- 
sation occurs when /i/ > m^r while kaon condensation is favored when fis > and niK < JTiyr- 
Using (j2ip and (j22p the condition (j37p can be rewritten as fj,j = m^r. For y = the bound- 
ary between the normal phase and the kaon condensed phase (see Eq. (jlO|) ) can be written as 



13 



10 
8 

X 

6 



4 
2 


1 2 3 4 5 

FIG. 3: Phase diagram of the Dirac spectrum m the fi^—x plane. The dashed curve is the Unc x ~ |/^^i^^/E. 

fj,s = 'mx/'^- The boundary between the pion condensed phase and the kaon condensed phase at 
y = (see Eq. (p9]) ) can be written as m/< = tjIt^. The reason is that since ^/ = 2/i5 = fi, the 
condensed phase with the hghtest boson wih dominate. 

IV. MEAN FIELD FREE ENERGY FOR n REPLICAS 

In this section we derive the mean field free energy of the partition function (jl7p . It is obtained 
from the free energy of the static chiral Lagrangian, which is still given by the static part of (j23p 
but now U € U{2n + 1). The baryon charge matrix and the mass matrix are now given by 

B„ = (5;.,odiag(ln,-ln,l), (45) 
M = diag(zl„,z*l„,m), (46) 

where is the n x n unit matrix. 




14 



The ansatz for the vacuum configurations is given by 



U - -^2n+l,p(/5) 



COS Ql 



sm ai 



cos Or, 



sm Or 



smai 



cos ai 



sm Or 



cosa„ 



(47) 



V ... Ij 

where R2n+i,piP) is a rotation by /? in the 2n + l,p-plane for some choice of p. 
In the normal phase all = so that U = 1 and 



-Ci^at = -T.{m + nz + nz* 



For a/c ^ 0, the minimum of >Cstat is either at /? = or at /? = 7r/2. For /? = we have 



(48) 



U 



cos ai 



sm ai 



smai 

cos On sin a„ 

cos ai 



sm 



cos an 
1 



(49) 



and find 



4^,^ = -mS - ^[2//2f2(1 - cos2 ak) + S(z + z*) cosofc]. 



(50) 



k=l 



The saddle point equation for each of the replicas is the same with nonzero solution given by 

^(z + z*) 

This is the pion condensed phase with free energy given by 



(51) 



*\2 



'-'Stat 



o o Y.^(z + z*) 



(52) 



8/i2F2 ■ 

For /? = vr/2, the ansatz (j47p for the vacuum configuration depends on -p. To excite mesons with 
a nonzero isospin charge we necessarily have that 1 < p < n. For k ^ p the saddle point equation 
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is again given by (j5ip and is the same as in the pion condensed case. For k = p the saddle point 
equation is given by ()35p found for n = 1. For arbitrary n we thus find that the free energy for the 
kaon condensed phase is given by 

£1. = -E.-2„m'F^-(»-1)5^^^-5^|^. (53) 
A. Phase Boundaries 

It is immediately clear that the phase boundaries between the normal phase and the pion 
condensed phase and between the pion condensed phase and the kaon condensed phase do not 
depend on n. However, the boundary between the kaon condensed phase and the normal phase 
does depend on n. Also for general n the kaon condensed region is split into a region with a spectral 
density that grows exponentially with the volume, and a region where it decreases exponentially 
with the volume. Because the spectral density in the pion condensed region does not depend 
exponentially on the volume for n replicas, we also have the relation 

This implies that the boundary of the oscillatory region with an exponentially large amplitude is 
given by the equality of the free energies of the pion condensed phase and the normal phase and 
therefore does not depend on n. 



V. CONCLUSIONS 



The Dirac spectrum for QCD with dynamical quarks at nonzero quark chemical potential can 
be split into three regions. A region with no eigenvalues, a region with a constant eigenvalue 
density equal to the quenched spectral density, and a region where the spectral density shows 
oscillations. This last region is further split into one region with an exponentially large and one 
region with an exponentially suppressed amplitude. These regions have been derived from the 
chiral Lagrangian for the generating functional of the eigenvalue density. The generating functional 
includes a conjugate quark and this makes it possible to reinterpret the chemical potential as a 
nonzero isospin and strangeness chemical potential. The three regions of the eigenvalue density 
correspond to a normal phase, a pion condensed phase and a kaon condensed phase of the generating 
functional, respectively. 
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These results were derived from the leading order chiral Lagrangian in the mean field limit 
which is valid for m, \z\ <^ F^/S and ^ F. However, we expect that the Dirac spectrum beyond 
this domain will behave in a similar fashion. 

The mean field limit is the lowest order term in the standard expansion of chiral perturbation 
theory 



33( 1 ■ The results presented here thus form the basis for an analysis that contains loop 
effects as well as higher order terms in the chiral Lagrangian. The inclusion of loop effects in 
chiral perturbation theory at nonzero chemical potential is well defined 



35|, 



36|. 



34l ] and can be extended 



to include the effects of a nonzero temperature as well 
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